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SUMMARY 


If- is shown that under' cer tai h ■ e ondit i on s . a two-dimen- 
sional ■ sub s oni c comjpre ssible flow around an airfoil profile - 
can-be derived from an inc ompr e s sibl e flo w * ar ound another ■ J 
profile.- The connection between these ' tw6 “ c on jugate flows 0 • 
is given by a simple conformal transformation of the resped- 
five holograph planes. . 

The transformation of a given incompressible flow into •• 
a compressible' flow around a slightly ' dist ort ed- pr of ile - r a-’ * 
duces to the integration of a linear ' part ial differential, 
equation in the physical plane of the incbmpressible flow." 
An approximate solution of ' this equation is indicated'. ■ Fur-' 
ther research is necessary in order to extend the ' applica-. 
bility of the method and in order to reduce the computational 
work involved in the- rigorous' solution to an acceptable min- 
imum . • t 


The transformation of an incompressible flow into a 
compressible one- can be carried -out c omplefc ely 1 and' in a' 
closed form under the assumption of the linearized .pressure- 
density- relation . The final formulas' represent an extension 
of the result of- Von Karman' and T si en , ■ t o ‘ whi ch‘ they reduce 
in the special- case of a flow without' circulation. .It. is. • 
shown ’ that essentially all compressible flows can- be obtained 
by this method. - ‘ „ 
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INTRODUCTION .v v 


The high level <whiqh has "be ©p, *' 3t taclli'ed by the theory of 
t wo-dimensi onal incompressible flows is due to the fact that 
this theory is based upon a highly developed mathematical 
theory, that o:f analytic functions of :A' complex variable. 
Fvery analytic function yields a- possible flow pattern and 
vice versa. - Furthermore, the use '.of ‘transformations per- 
formed by means of analytic functions (conformal transforma- 
tions) permits the derivation of all- possible flows from a 
few simple standard forms. It seems obvious that the theory 
of two-dimensional compressible flows (at least as far as 
subsonic flows are concerned) requires the development of a 
similar mathematical background. 

The theory of' si gma-monogeni c funct i on s (reference's 1 
and 2) is an. .attempt to study a class of. complex functions 
the role of which in gas dynamics is comparable to that of 
analyt.ic functions in the theory of incompressible flp.ws. 
Gelbart (reference 3) has outlined the application of. .this 
method to the study of ’ compressible flows. Further applica- 
tions depend upon the investigation of singularities of 
sigma-monogenic. functi ons . (Such an investigation is being 
conducted.) Reference also is made to a recent report by - 
Garrick and Kaplan (referenced) . The Investigation of 
transformations xvhich for the case of compressible flows take 
theplace of conformal transformations is the main theoreti- 
cal aim of the present report. • 

.'The following remarks may indicate in which way such 
transformations enter into the study of compressible flows 
around airfoil profiles. 

The differential equations governing the steady two- 
dimensional potential flow of a -compressible fluid are non- 
linear and therefore difficult- -to .treat as far a, s both the- 
oretical considerations _and numerical computations are con- 
cerned. .-Molenbr-oek (reference -.-5) and Tchaplygin (reference 
6) have shown that, the eonations become linear in the hodo- 
graph plane. There exist various methods of obtaining solu- 
tions of these hodograph equations, in particular of obtain- 
ing solutions which in a certain sense correspond to given 
solutions of the Oauchy-Ri emann equations - that is, to given 
incompressible flowB. Some of these methods are: separation 

of variables (so successfully used by Tchaplygin in solving 
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jet problems) , the method of .'integral operators (Bergman, 
reference 7), the method of sigma-monogenic functions (Bers 
and Celbart, reference l), and an approximate method of 
Temple and Tarwood , (r eference 8). furthermore, "by modifying 
the pre s sur e-density relation, the hodograph equations can 
he made to coincide with the 1 Cauchy-Bieinanh equations,’ as is 
done hy Tchaplygin (reference 6), Busemann (reference 9), 
Bemtchonko (reference 10), and, in a more rational way, hy 
Ton Karman and TBien (references 'll and: 12). 

HowB^e'r , the-real difficulty- lies not' only- in obtaining 
a solution hut in obtaining the "right” solution - that is, 
one which leads to a flow of a desired type in the physical 
plane, for example, to a flow' around a closed profile. This 
difficulty is illustrated hy the fact that even for the case 
when the hodograph equations' ar‘e ' simply the- ‘Can chy-Biemann 
equations, the computation of flows around closed profiles 
has until now heen carried out only for a special case (flows" 
without circulation). 

Therefore the study of ■ the ,- flow, in ■ the hp.d.ograph plajae 
must he supplemented hy the investigation of the mapping of 
the phy si pal . plane into -the', hodograph.. plane ‘and of- possible 
transformations of incompressible flows around closed pro- 
files into compressible . flow.S of the same type. ■ The -present ■* 
report is an attempt in this direction. 

.. The. methpd.s -outlined-. ih this, rep c5rt arp at present, re- 
stricted to flows which are everywhere subsonic. . Tlows of 
mixed type (subsonic main flow with locally supersonic re- 
gions) are .of more interest from the theoretical as “well as ... 
from the practical viewpoint. It is thought, however, that 
the solution of the problem of entirely - sub soni c flows is a 
necessary pr or equi si t e for .a sucoe s Sful .t be or et i cal .treatment 
of the much more difficult problem of mixed flows. 

This investigation, carried out at Brown University, 
was sponsored by and conducted with the . financial assistance 
of the National Advisory Committee for Aeronautics. The 
author i s -.in-deb t-ed t.o Mr : <J -. -B.- 'Diaz for • valuable assistance. 
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SYMBOLS 


a speed of sound 

ds ' line eleinent in the z-plkne; - line element of the 
profile P 

dS non~3uclidean length of ds 

da line element in the £~plane; line element of the 

profile n ■ ‘ * 

3(P) domain exterior to the profile P 

3 (f'l) domain exterior to the profile 

exp( ) exponential function of ( ) = e^ 

complex potential of an inc ompre s Bibl e flow in the 
I-plane, normalized so that G-'C 00 ) = 1 

holograph of a' compressible flow in the z-plane 

distorted holograph of a compressible flow in the 
z-plane 

imaginary unit; subscript referring to an incompress- 
ible flow 

imaginary part of' ( ) 

bound for the rati o of ■ maximal speed to stream speed 
in tho conjugate iiicompr e s sible flow 

Mach number 

stream Mach number 

modulus of the correspondence between two flows 

subscript referring to the state of fluid at rest 

a function f such that )f/R| — ^$>0 as H — R x 

being some specified limit 


o(i) 


H 

H* 


Im( ) 
K 


M 

Mcc 

n 

o 

o (R) 
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4 


1 


0(H) 


a function 
. 2—^ Hi , 


f such, that [ f/2| remains hounded as 
R x "being some specified limit 


P 

P 

q. 

q.* 

<u 


pressure 

profile in the z~plane 
speed 

distorted speed 

speed of an incompressible flow 


q. 

R 

2p 

E n 

Re ( 
T 

n* 

u* , 
x, 
X, 
z = 

Z = 

y 

i = 
0 


speed of. a f i ct it i ous 1 compr e s sible flow in the t-plane 
positive constant 

radius of curvature of the profile P 
radius of curvature of the profile H I 

) real part of ( ) 

coefficient of t he • symme tr i zed holograph equations 
v components of the velocity 

v* components of the distorted velocity 

y Cartesian coordinates in the z-plane 

Y Cartesian coordinates in the Z-plane 

x + -iy complex variable in the physical plane of the 
compressible flow 

X + iY auxiliary complex variable 

exponent in the polytropic relation 

% + IT) complex variable in the plane of the incompress- 
ible flow 

angle between the velocity vector of the compressible 
flow and the x-axis 


i 


J 
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0^ angle "between the velocity vector of the incompress- 
ible flow and the |-axis • ' 

X logarithm of the distorted speed 

value of the distorted speed at infinity for Y = -1 

t, T| Cartesian coordinates in the £-plane 

! f profile in the £-plane 

p density 

p density of the fictitious compressible . flow in the 

£ -plane 

cp velocity potential 

X angle between a line element ds and a streamline 

4^ stream function 

O complex potential of the conjugate incompressible flow 

complex potential of an incompressible flow around 
the circle [Z’t =* E 

co the point infinity; subscript referring to the state 

of the fluid at' infinity 

7 5 ~ complex conjugate of ( ) ’ 

]( ) } absolute value of ( ) 

The unit s are chosen so that j? 0 (st agnat i o n density) 
and a 0 ( 3pe ed of sound at a st agnat i on point ) are both 
equal t o unity . * 
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. : _ ANALYSIS 

I. GENERAL CONCEPTS 
fundamental Relations 


It will "be assumed that in a compressible fluid the 
pressure p is a given Increasing function of the density p. 
The "velocity of sound .a ..is given by - ■ 


- a s = 5E. 

dp 

If -the flow is irr ot at i onal , it follows from. Bernoulli 1 s 
e quat ion 



Po 


that the density is a given 'function of the speed q. Since 
the preceding equation can be written in the form 


q dq + a 3 ^2- = 0 

p 

the Mach number M" = q/a is given by the relation 


M 2 = _ il ^2. 

: ' P ^-q. 

■The units will be chosen so that 

a 0 ~ » Po' 3 -'-* (l) 

o referring to the state of the fluid at 
equivalent to the introduction of dimen si on- 
9./ a s * . p/.p o ; * •’ 


(the' subscript 
rest) . Thi sis 
less variable s 
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The pressure-density relation used in gas dynamics is 
of the form " ... 


P 


v 

— A + Bp 


( 2 ) 


This relation includes the case of an isothermal flow, where 


■ • Y = 1 

and that of an adi abat i c flow .with 


(3) 


1 < Y < 1.66 


(4) 


(The standard value of Y for air is 1.405.) The value 
Y = 2 .corresponds to the analogy "between a two-dimensional 
gas flow and a flow of’ water in an open channel; . (Cf., for 
instance, Von Karman , reference 12.) In the foregoing 'Cases 


A = 0, , B = p c 


The differential equations of a potential gas flow are con- 
siderably simplified by introducing the linearized -pressure- 
volume r elat i on with 

• Y = ~1 (5) 


and 


A = p 


+ a 3 p 

t» 03 1 00 


B 



P 


a 

oo 


where the subscript eo refers to the state of the fluid at 
infinity, and , p^ , .have been determined according to 

the actual pressure-density relation (with A = 0, B = p 0 , 

Y > l). Using this relation amounts to replacing the curve 
giving the actual pressure-density relation in the (l/p, p)- 
plane by its tangent at the point (l/p^, ) . The linearized 

pressure-density relation has been introduced by Von Karman 
and Tsien (references .11 and 1^) , (and formerly in a less gen- 
eral form byTehaplygin (reference 6), Demtchenko (reference 
1 0) , and Bu genian, (reference 9). 

The' relations between p, M, and' q. .obtained from . 

(2) depend only upon Y (and not upon A and B)-. . For Y=1 
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for y £ 1 


■ p - -s“4 a /S 

M = q. 

P = (i - -,.^/cv-i) 

■ n a = - ■, ■=■ " 



so that for V = — 1 


and 


vl + q® 

m 3 = — af 

- 'I + q^ 




( 6 ) 

(7) 

( 8 ) 

(9) 

( 10 ) 
( 11 ) 

( 12 ) 


In figure 1, p is ploth-ud 1 as' -a 1 function of q (for 
V s «i, l, 1.405) . - i . • . 

ffor V = —I , the flow, is always subsonic. Jor V ^ l , 
the flow is subsonic as' 1'ong- as 


a < l s 



■ 2 . , 

V + 1 


(13) 


It will he convenient to use, the distorted sneed q* 
(first introduced hy Busemann (reference 9)) 



( 14 ) 
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q * = 1 for M = 1 


Por 7=1 


q* = 


vW 


1 -*Vl - Q 2 . 


(15) 


(16) 


for 7 > 1 

1 - J 1 - M* 


q* = 


1 W 1 - M s 

and for 7 = -1 


1 + ycv-w 1 ) / ( 7 -f l)y/ 1 ~ M 3 


1 ~ 7(7- - 1 ) / ( 7 + 1)^1 ~ M s 


/7 + 1 

y 7^1 


* - 


(17) 


( 18 ) 


1 + 


y 1 + 


(in fig. 2. q and q* are plotted as functions of M for 
7 = 1.405.) ■■ .• 

Por an incompre s sitls fl'bw p is constant, M = 0, 
and q* = q. 

Equations of Motion 

The x and y components of the dimen si onl e ss velocity 
of a t wo-dimensi-onal potential steady gas flow, u and v, . 
satisfy the condition of irr otat i onality - ' • ' ■ • " 


8 u 

3 y 

and the continuity equation . 


B v 
B x 


= 0 


3 (pu) + ■ 3 ( pv) _• 0 
3 x By 


These equations imply the existence of a velocity potential 
qp(x,y) and of a stream function \J/(x,y) so that 
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u = ^22* v = -^32 
dx By 


(19) 


and 


pT = , |4 f .. P u = ft 

■ Bx By 


( 20 ) 


It follows tliat 


order equations 


, 9c 0 

_1 B\j/ 



3x 

p By 



% 



> 

Bep 

_1 Bjjr 



By 

p 3s 

wJ 


‘ either \{/ 

or ' <p 

lead" s 1- 'to 

A ( p |m\ + 

A ( p 

IlSlN 

= 0 

0x V Bx/ 

Sy 

ay/ 



( 21 ) 


( 22 ) 


JL (k + JL. (1 AA.\ = o 

3x \p 3s/ ... By Vp By/ 


(23) 


She fundamental equations ( 21 ) to (-23) are- of a purely kine- 
matic nature and hold independently of the. .equation of state. 
If the density p is considered as a given function of 
space (p = p(x,y)), the equations (22) are linear . and always 
of the elli~ptic type, no matter whether the flpw is subsonic 
or supersonic* .' •••' 

However, the important case is that in which the density 
is a given funct i on of nr essur e and therefore also a given 
function of the magnitude of the (dimensionless)' velocity' 

q =«/ u 3 + v 3 '. (Of. preceding sec.) In this- case equations 
(2l) to (23) are non-linear (more precisely: quasi-linear ) . 

The velocity distribution in .a given flow is uniquely 
determined by the boundary conditions and by the functional 
r elat i on 
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This remark justifies the use of the relation Y = -1 for 
subsonic flows of not too high Mach number. For, replacing 
the actual value of Y by -1 changes the function p = f(q) 
(and the differential equation for ;p), q < 1, only 
slightly. (Of. fig. 1.) From the equations given in the 
preceding section, it follows that the change of p is of 
the order of magnitude 


Y + 1 
8 


4 

M 


Molenbroek (reference 5) and T&haplygin (reference 6) 
shoved that linear equations can be obtained by considering 
cp and \j/ as functions of q and 0, where 0 is the 
angle between the velocity vector and the x-axis: 


0 = ( tan"" 1 ) — 
u 


(24) 


The equations 


take the form 


3cp 
d q 


8cp 

ae 


d\j/_ 

Sq 


id- m s ) 
pq ae 



(25) 


These equations can be brought to a symmetric form by replac- 
ing the independent variable q by 

X = log q* 


q* being the distorted speed. By virtue of (14), (25) can be 

written in the form 



^se. _ 
ax 


8\j/ 

ax 


% 

ae 


( 26 ) 
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with 

i 

T . A - M s 

P 

For the isothermal case (Y = l) .. 



(27) 


(23) 


For the polytropic case (Y > •' 1 ) 


I - V 2 , _1— Z _ 

; ( 

For the case of the linearized equation of state (Y = ~l) 

H — 1 (30) 

and equations (26) are Cauchy-R.i emann equations.-* -In figure 2, 
H is plotted-. as a function of the local Mach number M (for 
Y = 1.405). It' should 'te - noted' that- T is a known function 
of q and therefore also of q* 

The main advantage of the symmetric form (26) consists 
in the fact that the symmetric equations are invariant under 
conformal transformation of- the 6, X-plane: If new inde- 

pendent variables | and T) are introduced by setting 


1 - 


Y - 1 


% 


(Y+l) /( Y-l) 


I + in = f(q + ix) 

F being an analytic function, then 

9cp _ _ 3\j/ '■ 

H a-H 


3cp - . rp 3\j/ 
ST) “ 3 1 


( 31 ) 



14 


K ACA TN IT o . 969 


with. 


T = T<jq jjq* ( 4 , T) ) j- 

In fact, (26) can he written in the form 

2$. il . Ise. _ m /±k §1 + ^ M 

3| 36 Br| 80 “ \di 3X 3rj 3AJ 


^32. £L + 2 se. §n _ _ T f^L 51 + m £n\ 
84 sx 8<n ax -\84 ae • an aej 


Eliminating the derivatives' of 
Cauchy~Ri emann equations 


£ -and T) h y means of the 


3 £ __ 8r> 3 £ _ 3 ti 

30 = 3X’ 3X = “ 30 


yields -( 31 ) ; 

In particular, the distorted velocity may. he defined as 
follows: 


u*'- iv* = e*“ i6 = q*e H * i6 . 


(3 2) 


Then 


8cp 

3u* 


3cp _ 


3 v* 


.T ‘Jiit 

3v* 


i it 

au* 


(33) 


Since T is a given function of u* s + v* a = q* s , the sec** 
ond- order equations obtained from (33) are 
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where 


Acp - 2 


T ' 

(u* JSEL + 

3cp \ 

T 

\ 3u* 

3v*/ 

T * 

(XL* M. + 

y * 

T 

\ 3u* 

av*y 


0 


0 


A = (8 3 /du* 3 ) + 0 3 /Bv* 3 ). and' T ' 


dT/dq.* 3 


Distorted Hodograph of a Flow around a Profile 

Only the following types of flow will he considered in 
this paper. The flow covers the domain exterior to a pro- 
file P (domain E(P))’l At infinity the flow approaches a 
uniform flow in the positive x-direction, so that 


lim 
, x 3 +y 3 - 


u = 


a ~oo > 


0 , 


CO 


lim 

x 3 +y 3 - 


v = 0 


00 


(34 y 


The flow is everywhere subsonic (15 < l) . The profile P is 
a streamline of the flow; P is a piecewise analytic curve 
possessing at most two’ sharp edges. If there are sharp edges, 
the Kutta— Joukowski condition is satisfied. There are exactly 
two stagnation points, both situated on P. A uniform flow 
u = constant, v = 0 is excluded'. (For the sake of mathemat- 
ical discussion it is convenient to admit as ’’profiles 1 ' P 
curves which intersect themselves ' in a finite number of 
points. The exterior E(p) is then a partly multiply cov- 
ered Riemann surface.) 


Incompressible flows considered will be subject 
same restrictions, except that the edges need not be 


t o the 
sharp . 


The transformation 

u = u( x, y) , v = v( x, y ) (35) 

takes E(P) into a domain H of the (u', -v ) -plane ; H (the 
hodograph of the flow) is, in gener al, mult iply covered. The 
t r an s f o r mat i on 
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u* = SLL u = q* . cos 0, v* „= ill y = (j* ein fl (36) 

<L d 

takes E(P) Into a domain H* of the ( u* , -v* ) -plane . This 
domain will he denoted "as tire distorted hod-ograph of the flow 
( q* and -0 are the polar coordinates in the ( u* , -v* ) -plane ) . 
(Cf. fig. 3, (a), (b), (c).) 

It is known that in the case of an incompressible fluid 
the mapping of tine' flow' into it’s hodograph is conformal. It 
will be shown that the mapping of E(P) into H* can be con- 
sidered as conformal if angles and distances in E(P) are 
measured by. means of a oertain fiiemanrian metric generated by 
the flow. 


Transformations Conformal v/ith P.espect to a Given Plow 

Given a subsonic compressible flow covering a domain D 
of the (x,y)-plane. Let ds = (dx, dy) be an infinitesimal 
line element situated at a point, x-, y of D , M the value 

of the Mach number at this point, and x the angle between 

ds and the streamline passing through x,y. The non- 
Euclidean length of the line element ds shall be defined as 

. dS * dejl"- M s sin 3 x (37) 

- ' . ‘ t 

Let a be the angle between the 'line element ds and 
the x-iaxis-. . Sinc.e. X- - *(6 - a,) and dx' = ‘ds cos a, 

dy = ds sin a, it is easirly seen that (-37) can be written 

in the ,f orm ■ 

dS 3 = e dx s + 2f dx dy + g dy 2 (33) 

wher.e . 

e = l*- M 3 sin 3 6 

• . f ' = M S - sin, 0. ccs '6 >' (39) 

■ g = 1 - “M 3 " cos 5 '0- 


Thus (37) is a Hiemann metric. 
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The non-Euclidean ' angle A between two line elements 
ds = (dx, dy) , 6s = (fix, Sy) ■ ■ situated at the same point 
x,y is defined by ; 

. e dx 6x + f(dx 6y + dy 6x) + g dy 6y 

COS A = . : 

dS 6S 

6S being the non-Euclidean length of 6s. 

A transformation 

| = £(x,y), “H = Tl(x,y) (40) 

of D into a (simply or multiply covered) domain A of the 
(i,Tt)-pl ane will be called conformal with respect to the flow 
if it preserves the pense of 'rotation and takes each non- 
Euclidean angle A (in the' '( x , y) -plane ) into the Euclidean 
angle A. An equivalent definition is the following. The 
transformation (4.0) is conformal with respect to the flow if 
it preserves the sense of rotation and if the ratio 



where dS is the non-Euclidean length of a line element ds 
in the (x,y)— plane and ■' do - , t he . Euclidean length of its 
image in the (£, T] )— plane-, depends only upon the position (but 
not upon the direction) of ds. The symbol A is called the 
local factor of magnification . 

If D is mapped conformally with respect to the flow 
into A, and A is mapped conformally (in the ordinary sense 
into A* i the resulting mapping of ' E into A r is conformal 
with respect to flow. Conversely, if D is mapped conform- 
ally with respect to the flow into, both A and A‘» the re- 
sulting mapping' of A into' A 1 is conformal in the ordi- 
nary sense. 

Any trans format i on - (40)' pjf E(P) into the ( £ ,T1 )- plane 
i s conformal with r e Bpect t o the flow if the potential 
funct 1 on cp and the stream func 1 1 on \J/ considered a s func - 
tions of l and T| sat i sfy the differential equations ( 3 1 ) . 

This. follows from lemma 1 proved in the appendix by set- 
ting A = l/p, B = T . 
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If the. foregoing result is used, and it is noted that cp 
and satisfy equations (33), the ■ foil owi ng important theo- 

rem is seen to 1)6 true: 

The mapping , of the exterior ■ . E.(P) of a. profile P into 
the di st ort ed h-odogr aph of a subsonic compressible f 1 o w 
around P i_s c onf or mal with , respect t o this flow . 


Mapping of a Compressible Plow into a Domain 
Exterior t o an Arbitrary Profile 

■ i 

The distorted hodograph of a flow around an airfoil is 
a simply connected Eiemann surface bounded b'y a closed curve 
(the image of P), By a known theorem of funct'ion theory,- it 
is possible to map H* conformally into an arbitrary simply 
connected domain. -Therefore, it is possible to map S(p) 
conformally, with respect' to the flow into the domain. 3(H) of 
the ( t ,T]')-plane, .exterior to a given profile ~i . ' The mapping 
can be cho.s.en so that- the point % +111 = 00 corresponds to ' 
the point x + iy = oo and that the horizontal direction at 
infinity is preserved (i.e., at oo the direction parallel to 
the x-axis is taken into the direction parallel to the £~ 
axis). ^Furthermore, by eventually changing the size but not 
the shape of H it is possible to obtain- a mapping for which 
the local factor of magnification is equal t o 1 at infinity. 
Since at infinity the metric (37) approaches the metric with 
constant coefficients 

dS 3 = dx 3 + (1 - M^dy 3 


the above conditions mean that 


t s + ri 2 -fr-co, |£— >i, |£-»o 

ox 3y 


M'— > o, M. 

dx By 






M, 


CO 


(41) 


as 


x 2 + y s - 


00 


A transformation satisfying these conditions will be called 
normalized . 
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If E(P) is mapped conformally with respect to the flow 
into B(n) , the resulting correspondence "between the. .points 
of H* and those of E (l~0 is conformal. Therefore, 
u* - iv* is a one-valued analyt ic * functi on of | + iT) and 
can he developed in. a Laurent series for sufficiently large 
values of [t + iT) J . ■ The following result will he used later 

Lemma. 2t If the manning of E ( P ) i nt o E (H) i_s nor - 
mail zed , then 


,u* — iv* — a* — 

*■3- 00 


ir T ,i ■ 

2tt X + i'H 


where 


and 


r x 


m® sJk, r 

lco ' ' 


(42) 


(’43) 


r = j U dx +' V, dy 

1 s the cir culat i on o f the c omnr e s si hi e flow . 
The proof-. will he found in the appendix. 


(44) 


Conjugate Blows- 

» v 

G-iven a compressible- flow- (in: the (x , jjr) -pl.an.e ) around 
the profile. P. and an incompressible flow (in the . (|,T))- 
plane) around a profile" O'. The complex- potent ial of the in- 
compressible .flow will- he -deho-ted. by ft(£) = cp ± + ity, . Its 
complex velocity is" ... 

■ u i - i- v i = q.ie" t6 'i = Mkr- > ] 

■ . 

Since it is assumed that at infinity 6i ‘= 0* O can he 
written in the ‘form 

0(5) - 4 lfCB G(5). = 1 . * 

I 

The two flows will he called "c on .iugat e ( modulo n), if 
there exists a r.eal number . n, 
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‘ . . ' ' 0 ‘< a < 3 

such that' 'the transformation' 

. .. , ,u j - iVj_ = (u* - ‘iv*) n (45) 

takes the distorted hodograph H* of . the compressible flow 
into the hodograph of the incompre ssihle flow. The con- 

nection between conjugate flows is shown in figure 3. 

The mapping. (45 ) defines a mapping of - E(P) into E (fl) . 
This mapping is conformal with respect to the flow around P. 
Por so is the mapping of E(p) into H* ; and the mapping of 
H* into Hj_ (given by (45) as well as the mapping of Hj_ 

into E (fl) are conformal- in. the ordinary sense. The mapping 
of E(P) into E(m) preserves th/fe 'point infinity and the 
horizontal direction at infinity (for at infinity both flows 
are horizontal). On P (on H) &( 6^ ) is the slope of the 

profile. According to (45) the. slopes at corresponding 
points are connected ;by the-.r elation 

6i = n& (46) 

(The slope is defined as the angle between the tangent to the 
profile and the positive x-axis; the tangent pointing in the 
direction of the flow.) 

Conversely, if it is possible to map E(P) into the do- 
main ‘ 3 (l~i) i n- the, ( | , T| )-plane 1 ext eri or to a profile H, by 
a transformation' which is conforma-l with respect to the com- 
pressible flow around . P., which preserves the point co and 
th o' horizontal direction at infinity, and which changes the 
slope of P according to (46), then the flow around P ie 
conjugate (modulo n) to an incompressible flow around H 
which has stagnation podnt-s a,t the. points into which the 
stagnation points at P are taken and the direction at in- 
finity of which is horizontal (provided such a flow exists). 

\ ■ 

Eor, let G(£) be the complex potential of sjich. a flow, 

C 1 ( co) = 1 . Set 

u i “* iv i = CL*® a '(0 

The mapping of E(fl) into H* is conformal - that, is, 

u* - iv* is a one-valued analytic function defined in E (H) . 
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Therefore, 6 _j= -Im 10g(u* iv*) is harmonic in E(n)j n .0 

coincides- on 1 I • with the. harmonic function 

Bi = -Im log(uj_ - iv^ ) . Therefore, (46) holdB throughout- 

3(l 0 . Since n log q* is conjugate to n 0 and 

log -q_i = log l>i - ivi 1 to 6 i 

log q j. . = n log q* + constant 


The above constant must vanish, for at infinity 

• ■ 1*® n * 4i,o, 

Therefore : 

. q.* n = q.i (47) 

and by (46) and (47), (45) holds. 

- If an incompressible flow is conjugate to a subsonic 
compressible flow, then 


‘ ’ 11, -flax K . ■ 

1 i , co 

wher.-e. q^; max 1 & the maximum speed. .and K depends, upon 
and the stream Mach number ■ "Moo of the compressible flow. 
For Q.*max must-be less than 1 ( cf. (l5)‘) and therefore, 
(47) 

1j , max _ /I max N ^ ^*-.n 
a a n g * * 00 

,05 ^ CD 


(48) 


n . 
t>y 


(49) 


Note that 
Moo. 


q* is a function of ’ q_ 


and therefore also of 


If the Incompressible flow with the complex potential 
Q( £) is conjugate to a given compressible flow 1 , so is the 
flow with the complex potential 

X nUO* A > 0 
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for ft has the same hodograph. ' Thus the conjugate profile i • 
can he chosen so that the mapping of E(P) into E(H):is 
normal! zed’. . . 

If the incompr e s sihle flow around. H conjugate to.' a, ' 
compressible flow around P is known, then the velocity dis- 
tribution (and therefore the pressure distribution) around P 
can be immediately computed. Por (46) gives the correspond- 
ence between the points of l~l and P, and the speods at 
corresponding points are given by (47). Eowever, the present 
method has not been developed sufficiently to permit a solu- 
tion. of the direct problem; to find the incompressible flow 
conjugate to a compressible ' flow around a given profile P. 
The following sections contain the solution of the inverse 
problem: to find a compressible flow around a closed profile 

conjugate to a given incompressible flow, and the discussion 
of the existence of conjugate flows, which is by no means 
self-evident . 


II -■ CONSTRUCTION OP SUBSONIC PLOWS AROUND A PROPILE UNDER 

THE ASSUMPTION OP THE LINEARIZED EQUATION OP STATE 

Simplifications Resulting from the Assumption *V = -1 

Throughout this chapter the pressure-density relation is 
assumed to have the linearized form - that i's , *Y -is set 
equal to -1. Under this assumption it can be shown that under 
'very 'general condit i on s .each .compressible flow possesses a 
conjugate incompressible flow and vice versa. The inverse 
problem can be solved completely and in a closed form. 

The assumption Y = -1 implies the following simplifi- 
cat ions. 1 

1. The differential equation of the potential cp in the 
physical plane takes the form 

(;.+ v s ) guv lls B_ + (1 + u 3 ) = 0 

Bx 3 ,Bxay _ 3y 3 

Bx By 


u 


(50) 
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Cp^ •• (3J3) andj-.O-.p);. ) Thus- .the- surface 
* • • / , *• ■ . . ... - 

z - CfxCx.y). . { 51:) 

(potential siirfa'ce ) : i's ! a minimal surface (surface of vanish- 
ing mean curvature).' '' 

2. The n on-Euclidean .length dS of a line element ds 
in the (x,y)-plane (cf. sec. Transformations Conformal with 
Respect to a Given Blow) ‘becomes equal to. 

dS = di */ 1 - M 2 

where dl is the Euclidean' length of the line element on 

the potential surfaces the projection of which is ds.. . , 

* * * 

. . r Bor,, if. ..t.he, angle . ‘between the direction of th.e velocity 
.vector and ; ds is denoted- "by - x, it follows -from (37) that 


. dl^ = (*1 + u 3 -) ■ dx 3 +- 2uv dx dy + (l .+ v s ) dy 3 
. .= dS 3 + (u &x-+ v dy) 2 • 

a ds s + ds 3 q 3 cos 3 \ 

• ds 3 ( 1 - M 3 sin 3 xV ' ' 

1 M 3 . 

_ - - dg 3 -v -• •• 

i; '.. - i”- ‘ - 

j. . • 

3. The term cp + i\J/ is an analytic function of u* - iv* 
and of the pompie-x var-iahle in any plane into -which 3(P) is 
mapped conformally w.ith respect t <? the flow. (Of. (30) , 

(33).) 

Bxi sfcence of ' Con jugate Blows 

It will he shown presently that to any compressible flow 
around a profile P in the z-plane (obeying the linearized 
equation of state) there .exists a conjugate flow of v an , incom- 
pressible fluid around a profile in the £ -plane, provided 

either of the following two conditions is satisfied. 
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(a) The compr e s sibl e flew is circulation-free. In this 
case the modulus n can he chosen at random. In particular, 
it is convenient to set n = 1, 

..... O) The compressible flow possesses a circulation and 
the Mach number at infinity is restricted by 

M. aj </374 = 0.866 . . . (52) 

' * ‘ * ‘ 

In this case the modulus n i.s given by ■ ... . 


(Of. fig. 4.)' 


n 



(53) 


Pgr thji proof, map ' E(P) .into the exterior of a circle 
X +'■? ■= E in the Z-plane (Z = X + iY, R being a conven- 
iently chosen constant) by abnormalized transformation which 
is conformal with respect to the given flow. Then the corre- 
spondence. bet ween the di st or.t ed . hodograph H* of the com- 
pressible flow and the domain [ Z j > E, is conformal - that 
is, u* - iv* is an analytic function of Z. By lemma 2, 
this function has the form 


u* — iv* = q* - --. r A . 3l + . . , (54) 

. y -oo .jbtt.. Z 

being given by (43). Furthermore,-' cp + 1\|/ is an ana- 
lytic function of u* — 'iv* and therefore also of Z, 

[Z] > E; \|/ = 0 on 1 Z | = E for \|f = .0 on P. If Z goes 
once around the circle JZ] = R, cp + i\J/ increases by P. 
Next , 

lim ^SE = lim + 3cp B 

. Z Z — > coldx dX By dX J 00 

lim |3B = lim JiSB J* + Jcg Sy] = 0 
a 00 9 Y z — > co • Vax 3 Y ■ By BYJ 


for the mapping is assumed .to. be hormaliz'ed. It follows 
t'hat . ... • ! . / 
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cp + I\|/ = qj Z + 


(> 



iT log % + constant 
■' 2tt . 


a x (z) 


Thus cp + i\|/ Is in tlie Z— plane the complex potential of a 
flow around the circle LZ'j = R. This flow necessarily pos- 
sesses two stagnation points on the cirqle iZl'-'E? namely, 
the images of the stagnation points on P. IPor at these 
points the lines \J/ = constant intersect the circle. It 
follows that ...... 


&.(cp + MO = . A \ IT 1 

dZ lao \ “ Z 2 / 2tt Z 


(55) 


vanishes at the same points (Z = S 1# Z = .'S 2 , S x = ~S 2 ) as 
does u* - iv* . ; : 

Assume that there exists an • incompressible flow around 
a profile H in the £-plane which is conjugate (modulo n) 
to the compre ssihle flow around P. Let 

z = z( D 

map E(H) conformally into [ Z 1 > E taking . £ = o° into- 
Z = cd. Without- loss of generality it may he assumed that 

Z' ( co) = 1 : ; 


Then the complex potential £) *= + ii^ of the incom- 

pressible flow must, in the Z-plane , he of the 'form 


n = + i^) • 


The complex velocity -of: tte con jugate flow, uf 
is given hy ’• 

U 4 - iVi = ■ 


- ivi, 


u i ~ ^ v i = ( u * ~ ’iv*) 11 


On the other hand. 
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Therefore 


ii . 5Q /aa - 

dZ dZ/ d£ 



.fig ±-lii / ( u * 


dZ 


i v* ) n 


so that 

» • f - 

Vt(z) - 



R a /Z a - ir/giT q^Z 
(u* - iv*) n 


dZ 


’(56) 


The numerator of the integrand in the above formula 
possesses simple zeros at the two stagnation points. It can 
he easily shown tha-t : at the two stagnation points (in the Z- 
plane) u* - iv* ' vanishes of an order not higher than 1. 

For at these points Im log(u* — i v * ) = -0 possesrss jumps 
of magnitude e^Tf, £ j/rr and e jgir being the angles at the 

stagnation points on P. Therefore Re log (u* iv.*} = log q. 

behaves at these points as g^loglz ~ ) . And it has been 

a s s.umed.. that P possesses .only sharp edges if' any, so that 
0 <■€!■<■ 1-. It follows that at'. S r and S 3 the integrand 

becomes infinite of an order less’ than 1, provided 0 < n < 2 
At all other points Z, |Z) > R, the integrand is different 
from both 0 and co . Eence th'e integration can be performed 
also along the circle IZ] = R. 

Now, since ' * was assumed to be a closed curve, the 
integral (56) must vanish if the integration is performed 
along the closed circle fZJ - R. By (54) 


1 

(u* - iv*) n 



+• 


a ££i. 1. 

^ 2-rr Z 


+ 



so that the integrand in (56) equals 


, l. + JL\ 11+;. . 

2tt z 

In order that the integral 1 taken along a closed 
vanish it is necessary and sufficient that 



'CO 



: 0 . ■ ■ ; 


cxirve should 
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or, by virtue of (43 ) , that •• 


(nA - mC-. I)r = 0 


(57) 


This’ condition, will he satisfied’ if either F = 0 or n is 
determined according to (53). In. this last case the condi- 
tion n < 2 yields the hound (52) f.or lio,. 

, .,«« ,v ■' ‘ ' * * ' * 

Conversely, if (57) holds and n<-'2, the function (56) 
maps [Zj > R conformally into the exterior of a closed 
profile M in the £-plane. (iTote that the derivative of 
this function does not vanish.') Thi s- 'mapping satisfies the 
conditions 


..." & (»)' = “» . ’ ( £.) = V 

' e ‘ ... ‘ 

The resulting mapping of . E(P) into E(H) , is conformal 
with respect to. the flow .around.. P. Obviously, . . 


-• . • 

CO w _ 

, <t . r ...... v . • 

is the oompleX' potential of an incompressible flaw around ‘t l . 
The .complex velocity of this flow is equal to (u* - iv*) 31 . 
Thus the -flow . is conjugat-e to the compr e s.sible flow around P 
and the . assertions- formulated in- the beginning of this sec- 
tion are proved. • •_ • 


Properties of the Conjugate Plow 

Prom the construction of the conjugate flow given. in' the 
preceding section, it follows that the circulation of the 
conjugate incompressible fJLow is equal to;-... 



(r being the circulation of the compressible flow), provided 
n has been chosen so that the mapping of E(P) into E(H) 
is normalized. 
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For the conjugate incompressible flow- the rat'i'o 
9.* ma J<lA cannot exceed K = q*- 11 . (Of. end of pt . I.) 

X | utc& X 1 f CO CO 

lor Y si -1, 

- X ± J 1 - ^ 


K = 


-SP for n = 1 


M 


00 


V. 


K = 


r • 

+ Ji '/x-mIo 


< 5-8) 


M'm ■ 


for n = 1/Vi - M*, 

f -J 


(The yalue.s ■ of K' are plotted, in fig. 5.) 


It should be noted that the profile 1 I constructed in 
the preceding section is necessarily; closed but need not be 
simple - that is, i t might intersect itself, in which case 
E (n) would be. partly multiply covered. .For this reason such 
physically impossible flows were included in the di ecus si on.. 
It is easily been that H 'will always be simple if P is 
c onvex , 


Solution of the --Inver so Problem 

Suppose that an incompressible flow around a profile 1 1 
in the' £ -plane is- giv'en and it ‘ is known that this flow is 
con jugat e” (modulo n) to a compressible flow around a profile 
P ■ in the' z-plane . ■' Thi s section contains the derivation of 
th'e formulas which- permit finding P and the compressible 
flow around P. It will turn out that these- formulas always 
yield a compressible flow around a closed profile, even if 
l“l and the flow around H are chosen at random. 

Let the complex potential of the incompressible flow be 
given in the form - 

fi(£)'= q ijO0 G(£'), G-'Cco) =1 

Since cp + i\|/ (complex potential of the compressible 
flow^ considered as a function £ is' analytic and real 

on II and since d(rp + i\|/)/d£ vanishes at the stagnation 
points of the incompressible .-flofr, 


cp + i\|/ = 0G( £ ) 



* 
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0 "being a positive constant. 'Without loss of 
may he assumed, that, the mapping pf E(P) -into 
normalized. * Then ' •*' "" ‘ 

generality it 
SI (PI) i s 




8 

l! 

o 




*Tow, let ds he a line element on P and dcr the 
corresponding line element on . O , Then 




^32 - 0 d& = —Is 2_ Q . 
,, . A& Act q i(0 o 1 




On the other hand, on P 

a 


* 


dqp 

dl"-* 



• 

Since the increase of cp on ds is eq_ual to 
ep on d®", and since == q* 11 . 

the increase , of 


A 

By (18) 

d - s _ - ...a* . <ii - q* n 

4 ■ 4 ._ 



* 

» • i • r 

= i h • 

Q. 2 ■ 




so that 

> * » • * 

■ i 4 ; *» . 

*r» ‘ 



> 


ds q*n-i 

: : dcr ^r.^n+1 . .. . 

# 


* 

' ' ; . 
i . 

’ ,„i-x/n ■ - i+i/n 

_ *i ~ *i ' 

a *n-l _ a *n+l 

: , ■ , r 1 “ ^co • co • 

(59) 


'A 

Since 

’ B r • 

■? f t 

ai = q i>OT l&‘* c't ) l = 4* tt r&’ (O'! 

, 


* 

• 

• 

.(■59) can he 

•• ■ ’ * 

written as .. . 
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If = ( 60 ) 

where p, =c q* - that is, 


p, = f ig-- (61) 

1 + Vl - M|, 

Next, let dz and d£ he the complex line elements on P 
and n, respectively; 

dz = ds = ds e 3- 0 i / n. 

d£ = do" e i6i 

Then 

» : ►.* . 

iz = dff e 10 l/ n (63) 

do 

Since 

|G'(£)I = &' (£) e i6 i = G f ( £) e~ i6 i 


there is obtained by (60) and (63) 

dz = d£ - y, 3 &• ( 0 1+X//n ^ A 

' 1 - H 3 L J 

Integration yields the fallowing representation of the pro- 
file P: 

z ~ constant ^ ^ d £ “ H. 3 J' G-'(£) 1+1 ^ a d £ j* (63) 

the integration being performed along H, (The value of thQ 
constant factor affects only the size of P.) 

Por a circulation-free flow and n = 1 . this formula 
simplifies to '■ 

z = constant U ~ l^ 3 T &,S d£ \ (64) 
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(64) and (6l) are exactly the formulas given in Tsien's 
paper (reference 14). 


Parametric Representation of Subsonic Compressible Flows 

Formula (63) has been derived under the assumption that 
the integration is being performed along l l and that the 
existence of a conjugate compressible flow is known before- 
hand. Both .conditions are unessential. For the following 
general result holds: 

Let 0 ( £ ) be. the conrplex not ential of an incom-pr essiole 
flow ar ound a. vr o f i 1 e H | in. the £ — -plane . 

&' M = 1 (65) 

Let be a. real number- such that 


0 < <1 if G is one-valued, 

0 < M < J 3 / 4 if Or is multi-valued. 


Set 


> = 


K 


CO 




n 


Co 


[ 1+6, 0 < 5 < 1, if G is one-valued 


n = 


l//l - M cq if G is multi-valued 


( 66 ) 


( 67) 


-n 


I f f or £ e E (i i) , 


then the funct Ion 


E tr |J, 


G 1 ( £ ) 1 < E 


= nf r cifn 1 ~ 1 / n 


z = °| / *'(£) 7 H - / G»(£) A " i/n d£ J. (70) 


i + i /n 


} 


( 68 ) 

(69) 
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C > 0 , maxi s E (H) in a one -t o- one manner into the ext er i or. 
E ( P ) o_f a, closed xir of i le P in. the z-xla ne t aking £ = oo 
into z • — > co , 

cp = 2Q H Re G(£) 

c onsider ed as a. funct 1 on of x and y is. the not ent ial .of. 
a. suh s oni c a omxir e s s i hi e f 1 ow ar ound P (obeying the linear- 
ised equation of state”) of stream Mach number . 

The t wo flows ar e c on .tugat e ( modul o n ) . The mapping of 
E(P) into E (Tj) i_§. normali zed by choosing 0 = l/(l - P> a ). 

The proof of the mapping properties of the function (70) 
will be found in the appendix under 0. Since the mapping of 
E(i~ t) into E(P) is one-to-one, cp may be considered as a 
function of (x,y) . Equati ons ( 70) and (7l) may be rewritten 
in the form 


x = Re f(£), 7 = He g( £ ) ♦ cp = Re h(£) (72) 

where the analytic functions f, g, and h are given by 


til) = 0 

g(^) 

h( l) » 20 p, a 


(G-i 1 - 1 / 11 _ p, s g. i 1 + 1 / n ) 


= — i 0 J* (ff » + ^ a' 1+1 / n ) d i 


Since 

f ,a + g' 3 + h ,S = 0 

(72) is the well-known Tfeierstrass parameter representation 
of a minimal surface. Hence, cp(x.y) satisfies equation (50) 
and therefore is a potential of a compressible flow. 


It is shown in the appendix under C that as £ — > co 


dx 


~5> C(l 



dx 

Sr) 


~ > 0, 




—> o, 


^ 0(1 + V. 3 ) 

3r) 
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(Of. (03) f ) ?y (66) and ( 71 ) title implies that as z — > o° 


dcp 

9x 


-§tL 


-> - — > 0 


SSB. 

dy 


Thus, at infinity the compressible flow is parallel to the x- 
axi s and 


2p, 


VL 


^ 1 - M> S Jl ~ M 


a 

CO 


so that Moo is actually the stream Mach number. 


It is also shown in the appendix under C that the direc- 
tion normal to H is taken Into the direction normal to P. 
Since the normal derivative of cp on H vanishes, ,so does 
the normal derivative of cp on P. Thus P is a streamline 
of the compressible flow. 


The fact that the two flows are conjugate follows simply 
by comparing the velocities at corresponding points of P 
and f~l . The details of this computation may be omitted. 


Prom the preceding section it follows that the paramet- 
ric representation (70) and (71)' yields all flows satisfying 
the conditions stated in the beginning of part II. (Note 
that neither H 'nor P are necessarily simple curves.) 


Construction Of a Compressible Plow around a Profile 
Similar to a Given Profile 

Suppose a profile. ?! and a point S on this profile 
are given and it is desired to find a subsonic compressible 
flow around a profile P- similar to P 1( possessing a pre- 
scribed stream Maph number < Jz /4) and a . stagnation 

point near the jp^int S. (Since S is determined by the 
angle of attack file last requirement determines approximately 
the position of with respect to the undisturbed flow.) 

This problem qh?} ^e solved as follows. 

Por the sake of definiteness it will be assumed that the 
profile ?! has one sharp (trailing) edge. It may be as- 
sumed that the function mapping P x conformally Into a circle 
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is known. In fact, this function c&n be easily computed.- 
(See reference 13.) The first step consists in forming an 
incompre ssible flow around P x which possesses stagnation 
points at the trailing edge and at 3.' The direction of 
this flow at infinity is taken as the | -direct i on in the 
plane of P x (£~p.lane). Now, if E(P X ) is mapped into 
)3| >1 by an analytic function ; 

' Z = F(£.) (73) 

which satisfies the conditions 



Prom is determined the modulus n by (53). Now, 

let’ 0 X • and C 3 ■ be two circles' pas sing -through.the points 
e -ico • which ■ inter sect at ' the angle nrr » ' The 

function" ‘ ■ • . ‘ 1 - .*-■■■ . ; - 


' i •'» Ji(Z) 

* * , 1 * 

inverse to' (73) maps the infinite domain bounded by $n arc 
of 0 X and ' an arc of C s into the eStefior of some closed 
profile f~7 . Profile n possesses two singular points; the 
trailing edge which c oincide s ■ wi th the trailing edge of- P x 
and the point S. The angles there are nf3 and nir, {3 
being the angle at the- trailing edge , of; . , P x . 


Since’ a domain bounded by two- circular arcs can easily 
be mapped into- the exterior of a circle, it is easy to .com- 
pute incompressible, flows around’ !~l . ■ ■ 

From M w is determined q^ (by' ('ll)) and q 1 ^ (by. 
(18)). Let G-(t’) be the complex potential of an incompress- 
ible flow around T*l. which' satisfies the Kut t.a- JoukowBki 
condition at the trailing edge and the condition 


G- 1 (cd ) 


1 
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It is easily seen that this flow possesses a stagnation point 
at S. The same is true for the flow with the complex poten- 
tial 

ITow a compr essihle flow around a profile P is con- 
structed which is conjugate (modulo n) to the above incom- 
pressible flow; P is given by (70). The velocity potential 
is given by (71). The distorted speed q* of the compress- 
ible flow at a point z of P is equal to 


9 * = 



q^ being the speed of .the incompressible flow at the corre- 
sponding point of P! . Prom q* , q i s 'det ermined by (18) . 


Since the angle between the x-axis and the tangent to 
P at a point z is equal to l/n times the angle between 
the |-axis and the tangent to H at the corresponding point 
£, it is, seen that the profile P will possess only one 
Bharp edge, the angle there being P. 

The profile distortion (i.e., the difference between P x 
and P) is due (l) to the difference between P 2 and I - ! and 

( S) to the difference between fl and P. This distortio.n 
will be small if is not too large. For then u. 3 is 

small and n close to 1. (See table II; in fig. 4, n and 
H- 3 are plotted as functions of M^,.) Therefore, the circles 
C lt 0 2 are close to the unit circle and f~l close to Pj. . 
Secondly, 1— l/n will be small and therefore the first term 
in (70) will be close to £ while the second term will be 
small as compared to the first. 


Of course, 

other ways. 

If 

point at S 

i s 

possible to 

set 


it is possible to construct I - ! in many 
the flow around P 2 which has a stagnation 
circulation-free (i.e., if a = 0), it is 
n = 1. Then Pj_ coincides with n . 
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. ' ■ Alternative . Formulas 

It is useful to write the formulas transforming an in- 
compressible flow Into a compressible flow in 'a different 
form. The (dimensionless) speed at infinity. will be 

used as the parameter characterizing the flow. 

.from formulas (52), (61), and (11) it follows that 

n = y 1 + q. 3 (74) 

r V ^ CD 


n.3 = 3L.H-1 

■ n +■ 1 


(75) 


The function mapping the profile D into the profile P 
takes the form 


z - 



S'(O l + l/n it 


(76) 


The potential at the cdmpre gsible flow is given by 


Cp - 'Jn 3 - 1 G 


(77) 


(The arbitrary constant appearing in ( 70) and (71) has been 
chosen as i 1 / ( 1 ~.p» s ) = (n + l)/s).) Finally, the speed q 
of the compressible flow around P is given by 


s (fvD i l 6l ^>! 1/a 

i - (B-.}.) bi(g)[»/ n 

^n + l'.’ ■ ' 


( 78 ) 


This follows, immediately from formulas (47) and (18),. by- 
noting that in this case 


q i = u n |G> 1 
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Tsien's simpler formulas which are valid for the case 
of a circulation-free flow may he rewritten similarly. 


n + 1 * n - 2 
z = £ : 


&'(£) s a£ 


.q. 


'n - I' 




(?9) 


p is given hy the same formula as before (formula (77))* 


III. CONSTRUCT I ON'. OP SUBSONIC PLOWS AROUND A PROPILE 
UNDER THE ASSUMPTION OP THE ACTUAL EQUATION OP STATE 
Existence of Conjugate Plows 


In this chapter it Is assumed that the presBure-density 
relation is not of the linearized form hut a general one, 
say the polytropic relation (V > l). It will he shown that 
the construction of a compressible flow for a given conjugate 
incompressible flow reduces t ;0 the solution of a boundary 
value problem for a linear partial differential equation in 
the phvsi cal plane of the conjugate flow. 

The discussion of the exi st ence of a conjugate flow for 
a given compressible flow around P can' be carried out in 
the same way as previously.. The essential difference con- 
sists in the fact that . so far it has been impossible to char- 
acterize completely all compressible flows possessing conju- 
gate incompressible f lows . and to determine a priori, the 
modulus n. 

Consider a compressible flow around a profile P in the 
z-plane . It has been shown that.it is possible to map B (Pi) 
into the exterior of a circle )Z| = R by a transformation 
which is conformal with respect to the flow’ and normalized 
at infinity. The two stagnation points of the flow are taken 
into two points 
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.Z v = Ee“ ia j Z 3 = Ee ~i(a+S) 


The flow will be. said to satisfy condit ion A if 

6 = tt - "2a 

that is, if 

Z = -Z 

S 1 


Assume that condition A is satisfied, and set 


r = -4 tt q^E. sin a 


(80) 


If 


0 < T /r < 2,/l - M 3 (81) 

Q CO 

(r being the circulation of the compressible flow), the 
flow will be said to satisfy condit i on B_, If r = 0, condi- 
tion B implies that P g = O'. 

In the case *Y = -1 , cp + i\[f 1 s an analytic function of 
Z; condition A is always satisfied and r is equal to r . 
Therefore condition 3 is satisfied for " 

T = 0 and for r f Q. M co < -sA 3 / 4 

It is con .lectur ed (but has not been proved) that A is 
always- satisfied and B is satisfied if does not exceed 

some limiting value. 

Condition A is' certainly satisfied (by reasons of sym- 
metry)- if either of the two geometrical conditions hold*: 

(i) P is symmetric with respect to the x-axi b and the 
flow is circulation-free. 

(ii) Both P and the flow around P are symmetric with 
respect to the y-axi s . (ThiB type of flow includes flows 
with circulation.) 

Plows of type (i) obviously satisfy condition B. In 
general, it may be assumed that B will be satisfied in many 
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cases. For the flow satisfying the linearized equation of 
state is a rather good approximation to the flow satisfying 
the actual equation of state and therefore r s should not 
differ too much from r . 

• 

I f a. compressible f 1 ow sat i sf i e s hoth conditions A 
and B, then it i s con.lugat e ( modulo n) t_o aa incompr e ssihle 
flow ar ound a profile n in. the t - plane . I f p = 0, n i s 
arbitrary and may be chosen as 1 . If. T ^ 0 , n i_s given by 



(. 82 ) 


Th.e proof is' almost exactly the same as the one given 
for the special case V, » ~1 and it wall suffice to sketch 
the argument . 


Let Qj.. be -the complex potential of a flow around 
! Z [ = R. which has stagnation points at Z-± and Z 2 and the ; 
velocity q^ at infinity,. Then Q has the form 


Q , = q (z + - i Ls . log Z + constant 

1 ■ 00 \ - Z/ 2tt 


r a being given by (80). As before, u* iv* is an analytic 
function of Z‘ and has' (according to .lemma 2) the form (42). 
If there exists a conjugate (modulo* n) incompressible flow in 
the " ^plane (around a' profile PI ), its complex- potential' fJ 
considered as a function of Z must have the form 

£2 = (X a real constant) 


As before, the function which 
be given by 

• -Z 



maps 


(u.* - 


I Z | > R into 


iv*)** 11 dZ 


E n> 


will 

(83) 


and the requirement that f~l be a closed profile leads to the 
condition 


nr i r s 
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or 



(84) 


which is satisfied if either 

• r = r - o 

a 

or a' is det ermined , by (82.). 


Conversely, if (84) holds, formula (83 
profile n (note that by.(8l), 0 < n < 2) 

preBsible flow around PI which .has the vel 

at infinity and stagnation points at £(Zj) 
conjugate (modulo n) to the given compressi 

It can be easily seen that the conditi 
only sufficient but also necessary for the 
conjugate flow. 


) yields a closed 


and t. 

he 

inoom- 


ocity 


t 03 ^ 

*n 

00 

and 


Z 3 ). 

is 

ble f 1 

ow . 



on s A 

, B 

ar e 

not 

exi st e: 

nee 

of a 



Solution of the Inverse Problem 


Suppose that an incompressible flow around the profile 
n in the £~plane is conjugate (modulo n) to a compressible 
flow around the profile P in the z-plane. It has been 
shown that <p and- a|/ (potential and stream function of the 
compressible flow) satisfy the equations (33) in the distorted 
hodograph plane. Since the correspondence between the dis- 
torted hodograph E* and E(D) is conformal, it follows 
that cp and \J/ , ■ considered as functions of £ and 7], sat- 
isfy the equations 


_ rn M 

3£ '3r| 

£& _ m 

drj , 3 £ 


(85) 


and the second-order i 

Jl J 

equations 

Ci 3cn\ JL J 

fl_ dgp~\ 

H 

It a£j 

r ’ ar^ 1 

it arj 

4j 

iil 



a£ 1 

L d ^J 

L aT iJ 


0 

0 


( 86 ) 


(87) 
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If the incompressible flow around f~l is known, T is 
a known function of £ and- T) . For. T 'is a known function 
of £ (see (27)) and therefore also of q* (see (14). And 
at corresponding points of E(P) and B ((~l) , q* and q^ 

(the speed of the ' inc ompr e ssible . flow) . are connected by the 
relati on 


*1 - 


(S8> 


Therefore the equations (85.) to (8?) are linear . 

The boundary conditions for and ■ are 

\l/ = constant, ™ = o 

d. indicating differentiation in the direction normal to FI , 
du - - 

At infinity cp and ^ must satisfy the conditions 


(89) 


and 


> C v > 0.’ 

at 1 a-n 


2i_>o. S4 
a I a-n 


> o 


(90) 


(91) 


where C x and G 3 are positive constants. This can be eas- 
ily verified by noting that the mapping of B(P) into B (PI) 
prese'r’ves the horizontal direction at infinity. The numeri- 
cal values of these constants are sof no consequence , since 
both the differential equations and .the boundary conditions 
are linear and homogeneous. 

.■ Function \J/ must always be one-valued; cp is one- 
valued only if T = 0. Moreover, cp and \p must satisfy 
the conditions: 


grad 3 ' cp < to. 'grad 3 < e» (9 2) 

Thus, it is seen that the equations (85) can be inter- 
preted physically as the equations of motion of a compressi- 
ble fluid of variable density p' vrhich';is a given function 
of space: 



42 


NACA TN No. 969 


P (I .Tl) = (93) 

J 1 - M3 

("I is a streamline of this flow and at infinity the flow has 
the positive ^-direction. 

Assume that cp is known as a function of | and T) 
and set 


q. = 




(94) 


(q is the speed of the fictitious compressible flow of den- 
sity p). Let ds he a line element on P and da the 
corresponding line element on f~l , Then 




so that 

d_s _ 
da “ q 


Since the angle between ds and the x-aiis is 0 and that 
between da and' the ^-aj-is'is 0^ = n0 , it is seen that a 
representation of the profile P can be obtained by 



(95) 


The connection between the profiles 
be expressed by the for mu las 


P and PI can also 


R 


R 


* _ n 5 

= n — , 

n * 


1 

a 


( 96 ) 
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where Hp - is the radius of curvature of the profile P at 
some po'int z and Hpj the radius of curvature of PI at 
the corresponding point £ . Angle 9 is known along n and 
so is q, for q* = / n . .If . q is known along Pi , then 

the profile P can he constructed graphically, using (95) or 
(96). 

Solution of. the. Inverse Problem (Continued) 

It has been shown in the preceding section that the con- 
struction of a compressible flow conjugate to, a given incom- 
pressible flow depends upon the solution of a classical 
boundary value problem for a linear partial differential 
equation (equation ( 36) )' "or equat i on • ( 8? ) . The integration 
may be performed not in E(n) but in a simpler domain, say 
in the domain exterior to the unit circle. For a conformal 
tran'sf ormat I Oh of • E(P!)' ’ into such a domain takes equation's 
(86), '( 87) ' int o equati ons of the same form and does not af- 
fect the auxiliary conditions. .Nevertheless, the actual in- 
tegration would be extremely laborious, especially in view of 
the fact that 'the coefficient T would be 'given -either numer- 
ically -or .by a very complicated analytical expression. Fur- 
• ther research is necessary in order to reduce the computation- 
al ‘work to an acceptable minimum. 

The physical interpretation of equations (85) given in 
the foregoing shows that these equations can be solved mechan- 
ically by G-. I. Taylor's well-known 'method of the " elect roi*. ■ 
lytic bath"., (See reference 14.) It should be noted that, 
whereas Taylor applied his method In order to obtain a se- 
quence of successive approximations to the solution of the 
direct problem (l.e., to the computation of the compressible 
flow past, a given prpfile), in this case the method immedi- 
ately fur ni she s. the exact solution of the ihverse problem. 

For slow flows (i.e., for flows where the local Mach 
number is small) the following approximate method can be 
used. 

, The .coefficient . T . in the equations (85) is equal to 1 
for M = p: and decreases very slowly as M increases to 
about 0..6, as seen from table I, where the polytropi.c reia- 
tionlwi.th Y = 1.405 has been assumed (cf. also, fig.' 23) . 

Therefore, for low Mach numbers, the density p "of the 
fictitious flow in the £-plane (see preceding section) is 
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almost constant so that q^ vill be a good approximation to 
q. Replacing q by q^ = q* 11 in the formulas (95), (96) 

yields 


and 



(97) 


Rp q*n g _ 1 

Rp 4 ’ 6 i n 


(98) 


Since is known as a -function of q^ and therefore 

q** 

known along f~l , the profile P can be immediately con- 
struct ed . 


Prom figure 2, where q and q,* are plotted as func- 
tions of M, it is seen that the profile distortion will be t 
small for small values of and for n close to 1. 

This approximate .method is based upon setting^ T = 1. 
The same assumption is made in the ? chaplygin-Karman-T sien 
approximate method. However, there the equation of state is 
changed accordingly and the resulting differential equation 
is integrated rigorously; whereas here the rigorous equation 
is solved approximately. (However, see also reference 12, 
p . 348 . ) 


Construction of a Plow around a Profile 
Similar to a G-iven Profile 

Suppose it is desired to construct a compressible flow 
around a profile P similar to a given profile P x and hav- 
ing a prescribed stagnation point S. For the sake of defi- 
niteness it is assumed that P x possesses one sharp trailing 
edge . 


Since the connection between n and is not known 

a priori it is advisable to star.t by choosing a value for n 
and constructing a profile f~l' such that if P is obtained 
from n, by using. (95) with this value of n a profile pos- 
sessing only one sharp edge will result. Profile FI can be 
constructed by the method described in the preceding chapter. 
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Let <G(£) be- the potential .of a flow around f~l which 
has the- Velocity 1 at infinity... It is now necessary to find 
a value q*^ such that if q* is determined as 



and T is determined as T‘ = T(q*), the values of q ob- 
tained bjr ' integrating equations (85) to 187) will lead to a 
closed profile P, when substituted in the relation (95). 

At the present state of the theory this can be achieved only 
by a tr ial-ahd-err or method. it is convenient to start with 
the value of q*, ? given by (74) and then change this value- so 
as 'to obtain a closed profile' P. Since this involves the 
integration of the equations (85) to (87) for different func- 
tions T, the" amount ' of computational work is rather consid- 
erable . 

*r 

In the case of' a circulation-free flow the situation is 
much simpler for q*, is Independent of n and n may be 
taken as 1. 

CONCLUDING REMAHKS * ■ ' - 


A method is given to transform a two-dimensional incom- 
pressible flow around a closed profile into a subsonic com- ■ - 
pr e ssib.le flow around another closed profile. The profile 
distortion is. small for small values of the' stream Mach num- 
ber . . 

In the case of the actual equation of .state the trans- 
formation depends upon the solution of a classical boundary 
value problem for a linear partial differential equation. An 
approximate method of solving 'this problem is indicated. 
Further theoretical work is required in order to establish 
the validity of the method in all cases and in order to reduce 
the amount of computational work. 

It 1 s • believed -that, after the solution of this inverse 
problem is completed, a way of solving the direct problem 
(computation of the flow around a given profile) will b,e open. 

If Ton Karmau'a and Tsien's linearized equation 'of state 
is assumed, this transformation is carried out completely and 



46 


NACA TN No 


9 69 


in a closed form not only for flows without circulation v 
(which wag already done by Tsien) hut for flows with circula- 
tion as well , 

Concerning the applications of the linearized equation 
of state ( Y = ~l) the following may be said. This equation 
of state can he applied consistently. But, it also seems 
worthwhile to try to use the assumption V = -1 only in 
order to obtain the values of the dimensionless speed q 'and. to 
compute the resulting Mach number by means of the rigorous 
equation of state. 

Other applications of solutions based upon setting 
y = ~1 also suggest, themselves, for instahee, to the solu- 
tion of the exact equation of motion by successive approxima- 
tions (by using tb.6 solution for. Y = -1, instead of the 
solution for an incompressible flow* as the first approxima- 
tion)-. 

Numerical examples and a comparison with other methods 
will' be 'given ’'in' a subsequent report* 

r 

Brown University, 

Providence, B. I., .April 29,- 1944;, ; 
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. APPENDIX 


A. This s.eption. contains, the proof, of the following 
Lemma 1 ; Let . ' 

i ’ 

t ~ i(x ,y), T| = T](x,y) (Al) 

he a traasf ox-mat , ! on of a domain D in the ( x. y) - plane into 
‘a, domain A .in the ( | , T) ) -- plane , Let 'qo and i{/ he. funct 1 ons 
of x and y defined in D . ( virtue of the above tran s- 

f ormat ion, they jn'ay also he o oh sidered as functions of | and 
I) defined in . A~T I_f and ^ sat 1 s fy in D the differ - 

ential equations 


- 

dx 


= A 


By 


•dgp''l 

By 


= -A 


ilk 

dx 


> A = A ( x, y) > 0 


( A2) 


a nd in A the differential eauat 1 o’n s - 


Be p. ■■ ■ : 




at B ,aW . 


>. B = BCe.n) > o 


(A3) 


Bc£' _ 

. B*n ” B a.t 


t hen , the tran-s format 1 on' (Al) jjj.' o onformal with r e spect t o the 
follow! ng Riemann metric defined in D; 


dS 3 * e dx s + 3f dx dy + g dy s 


(A4) 


where 
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e = 


f = 


cos 3 8 + sin. 2 0 

A 3 


"s 


C 1 - fi) 8in 


8 cos 8 




3' 


g = sin. 3 6 + — - cos 3 0 


(A5) 


a-od 6 is. the angle bet ween the line \jr = constant and the 
x~ axi s : 

tan 0 = M /9£ 

By/ Bx 


Geometrical -proof, (Cf. fig. 7.) - Equations (A2) show 
that the lines cp = constant and .\{f = .constant subdivide 
the (x,y)~plane into infinitesimal rectangles pf side ratio A. 
Similarly, equations (A3) express the fact that the lines 
cp = constant and \J / = constant subdivide the (| , Tl)-plane 
into infinitesimal rectangles o.f side ratio B. Therefore, in 
the neighborhood of some point (x,y) the mapping (Al) can be 
described as the product of a similarity transformation and a 
transformation which contracts all lengths in the direction 
of the line cp = constant in the rauio 3/A. But a mapping 
conformal with respect to t r he. metric (A4) is exactly such a 
mapping. Tor, let ds = (dx, dy) be a line element in the 
(x.y)-plane and let a be the angle between ibs direction 
and the x-direction. Then dx = ds cos a, dy = ds, Bin a, 
and x = * (0 - a) is the angle between this line .element' and 
the line = constant. A short computation shows that the 
non-Euclidean length of dB, as given by (A4) , is equal to 

dS = d.A - tl - B 2 /A s ) sin 2 x 


Thus, for ds parallel to the line cp = constant, 

dS = (B/A)ds, and for ds parallel to the line = constant, 

dS = ds . 

Analytical -proof .- Equations (A3) can be written, in the 


form 
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§£ i* + = B ££ + &l£ By\ 

Bx B£ . By B£ \Bx Sr) Sy Sr)/ 


Bcp Sx t Sep dy 
Sx Sr) Sy Sr) 

Sy use of the relations 


Vsx ag ay ag ) 


_ j 3 £_ j si. 

B| - J Sy' dr) “ Sy 


By _ 

H 


•J §a, Sx 

9 x Sr) 


J il 
Sx 


J _ B (x, y ) 

3 Cl , T) ) 

this can he written in the form 


an .. ase an = _ B ('14 ii _ ai M-) 

Bx dy Sy Sx ■• N 0 X ay ay ax/ 

dep d£ Bcg_ BJ_ p / B\j/ Br) . Syj/' Sn \ 

Bx Sy Sy Sx \ ax Sy Sy sx/ 


Introducing the values of S^r/Sx and d\J//Sy given by (A2) 
and using (A6) gives 


B 

-r COS 

6 il 

+ 1 

sin 8 

ii * 

sin 8 

in 

A 

Sx 

A 


By 

Bx 

sin 6 

B£ 

■ ■ 71. «M 

cos 

e ii 

3 

= — cos 

e in 

+ l 


Sx 


dy 

A - 

Sx 

A 


cos © in 
By 


Sy 
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Solving for B|/3x 



and 


sin 0 


B 


d£/3 y' gives 


cos 



0 .+ cos 3 



cos® 



(— 



sin 0 


cos 



If the notations given by (A5) are used, these equations can 
te written in the form 


it = 

3x 





f 3 



it = 

dy 


s 


in _ f in 

dx 5y 



But these are the well-known Beltrami equations which express 
the fact that the mapping (Al) is conformal with respect to 
the metr i c (A4) . 

R emar k ; The converse of the lemma is also true and can 
he proved similarly. A transf o'rmati on (Al) conformal with 
respect to the metric (A4) , (A5) takes cp and satisfying 

(A 2^ into functions which satisfy equations (A3) in the 
( I > ^1 ) -plane . 

B This section contains the proof of lemma 2. _I_f £(P) 
i s mapped into . £ (fl) by. a transf or mat ion 

I = £( x , y ) , T1 = Tl(x.y) 

which i s conformal with respect t o & subsonic compressible 
flow around P and such that 



5S 
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M 

bx 


1 . 


li 


e* + ti 2 


o, 


Si 


— 0 * 


SB 

Sy 


y 


1 - M 


CO 


(31) 


as x 3 + y 3 > as 


then u* - iv*, considered as a funct i on o f £ = | + iT| 

' has at the neighborhood of . £ = e» the ■ L-aur ent development 


where 


u* - iv* = q* - i + . . . 

• . 00 S'iT i 

r x -/i p 


(B 3 ) 


(B3) 


r being the circulati on o r f ■ the flow around P . 

The above transformation- can be written in the form 


i = A + [* 1 1> d x 1 

Bx ’ 


a 


+ / dy* 

ay 1 

V 


u To 


= A + ( x - a) + 


r Mill— >21 „ lldx 1 + •/ — 2 dy 1 

J 1 d3C ‘ i J* &y' 


a 

x y 


s rj = b d X . ,^aiK> . y, L > dy , 

a b 

x y 


= -B ; + ■/ 3M*I > b) > dx» + tf^y- b> + P (sMx^O _ N J dy , 




Bx 1 


,J b- 


v. . .3y 


where (a,b) is some point of ’E(P),' 


A = i (a.b) , ' B = Tl(a,b) 
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and 


I 3 

Co 


If use is made of l'Hopital's rule and (Bl), it is seen that 


I = 


■n = 



00 


Introducing 


shows that 


e 3 = x 3 + y s 


| = x + o(R) 

T) = ^ ^ y + o(B) 


(B4) 


Therefore 


1 1 


£ x + ill aa- y + o(R) x + i N q, y 


•© 


F “ °(^) ; 11 > 1 


Since u* - iv* i's a one-valued analytic function of £ 
in $(n) it possesses the Laurent development 


u* - iv* = q* + + a - s - + + . . . 

t t* 


(B5) 


■ M IQ 

Hence, for sufficiently large values of x + y , 


a + i£ 

u* - iv* = q* + — + o 

x x + iHocy NR 


0 


(Be) 
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so that 


u* = 5*,+ + o(l 


+ N 3 y 3 

r-n » 




_ aN y - px 


* +^0,7 


3 3 + °G) 


q* s = u k< 3 + v *' 


+ 2,. a ?. + . gM °> T 


+ 0 


* + ^a> y 


(l) 


a* = q* 


+ a. * ,.t Jlgg + „/.l 


X s + N 2 JT 
oo * 


<s) 


Next, q/q* is an analytic function of q* , so 


~t = + A <»(q* - ) + 0(.Jq* - q*} 3 ) 

* 4 00 

lq* - q* I 3 — > 0 


wher e 


since 


A = 

CO 


dq* N q 
1 


©! 


q * = c * 


# 


d,* 


<U* 


a * 

* 

j * q 


* = n* 


-•-*% Q- - 0 

q* ^co 7 


Af 

dq 


dq 


exp 


Vl - M ; 


dq 


- Sijl - H‘ 


(Of. the definition of the conjugate speed q* ♦ ) 


(B7) 

(BQ) 

(39) 

(BIO) 

that 

( Bll ) 

(B12) 
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By (Bll) and (BIO) .. ' *' 


A. = + A ax4 ~ $ll a y ‘ + -fl 

00 x 3 + N 2 y 2 
oo 


<1* <l£ 


Ther ef or e 


u = 


n* = a + -S_ r + 0 (T) 

" 4* »„ X s + Jfg y s W 


q 

v = - 2 *- V s 

q* 


= _sl y ,.z $*■■+ 

q* x s + N « y s \ R / 

00 


The 'ayttfcol p is an analytic function of q and 
fore also of q* , Hence 


P - P„ + V 1 * ~ 4^, > + °< I** ' " 

Id* - '11*1—8. O’" 


wher e 


B = -*£_ 

. ... “ r -dq* 

' : : ■ : V* = J <i» 

since, ,h;y :B:er n-oul 1 i 1 s equation, 


"..it*- 


By (B15) and (BIO) 


. £.k» 

dd d 


p = p + B + 

y Kco 05 .3 ... .«» „3 


<=> 


X" ■:+: Y* ■ ■ - >R ' 


(BIS) 

(B14) 

there- 

(B16) 

(B16) 

( B17 ) 


so that, hy (BIS), (B14), and (B17), 
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pu 


Pco «L» + 



ax + pN y 


K 7 S 


+ 



pv 


aN^ ~ Ex + 
q*' x 3 + N* y 3 ' 

CO CO " 



(B18) 

(B19) 


Nov, let 0 "be any simple closed curve containing the 
profile P in its interior. Then 


r = 


u dx + v dy 


(B20) 


0 


pv dx - pu dy 


c 


(B31) 


In particular, it is possible to take for C the ellipse Cj 
with the semi-axis 


the equation of which 


a = R, h = 
i s 

x S + S !L 7 s 

CO 


E 


N 


CO 


= p 2 


Por sufficiently large values of R the developments 
for u, v, pu, pv previously obtained may he introduced 
under the integral signs in (B30), (B2l). Then, using Green's 

theorem and denoting the interior of Cg by E^ results in 


r * 


%> 


dx + 


R" 



JL ( a3C + pN_ y) dx 


C R -Or 


’Ol*, 


■% (aNoo y - Ex) dy") + o(l) 


= - JL Bss. E JJ dx dy + o(l) 


* 


SP q 
N 


— ® TT + O ( 1 ) 
OO q CD 
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0 = - . j. .p q ' dy + -i- P -fp. iks (aN m y - Px)dx 

••• 1 “ " »* / l P » <£*. • .. 

^R ^R 


OO 


" Pj» 


Sje Rqq ( ax ’ *+ -y)dy\ + 6(1). 

J ' 



6.x dy + o(l) 


= " 2ap las tt 

q* 


+ o(l) 



Letting R 


whence 


oo yields 


r = . 5nL 3a, 


2rra,p c o 


isa 


P 


g«, r a*. 

2tt 4o» 



(BS2) 


If this is substituted into (B5), it can he seen that (B2), 
(B3) is verified. 

0, This section is devo’ted to the proof of the properties 
of the mapping function (70) whioh have been announced and 
used in deriving the parametric representation of subsonic 
flows with Y = -1 , 

First of all, this function maps f~l into a closed pro- 
file P and is one-valued in E(n) . In fact, -let J be a 
closed curve around PI and d > the increase in z . as £ 
goes once around J (J may coincide with (~| ). For suffi- 
ciently large values of ]£{, '$*(£) has the Laurent develop- 

ment 

G-‘(£) = 1 + + . . ., A real 

i£ £ 2 1 
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Therefore 



If Aj_ = 0, d = 0 . If A x 0, G(£) is not one-valued, 
so that n = ijt/'l - and therefore 

(l - l/n ) - n 3 (l + 1/n) = 0 

so that d = 0. 


<U = 

On n 


3 ft 

Next, let dz = dse 10 "be the complex line element on 
d cr the corresponding complex line element on f~l . 




O' ( O a \G'[ e” i 6i 


so that by (70) 


dz = Clff'l. 1 ^ (1 - 2//n ) e i9 i/ n d<T 


(Cl) 


( 02 ) 


and therefore 


ds „ i~i/u . s 3/u s 

— = C |G«[ ’ (l - p, ]G‘[ ' ) 

do - 


By virtue of (69) ds/du cannot vanish except at the two 
stagnation points where G' = 0, and. is positive elsewhere. 

Hence the . mapping of -fl int o B is one-to-one . 

Furthermore, it Is easily seen that z is finite for 
all finite values of £ and z = co for £ = oo. 

Finally, .the Jacobian . /. _ 

• d(y, y ) 

. a( i, n) . 

does not vani sh in E (fl) . r F or 
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(03) 


The expression on the right-hand side is not zero, for it is 
always assumed that there are no stagnation points within the 

flow and (l ~ p, 4 JG- 1 { -4f / n ) >0 hy virtue of (69). 

Prom the preceding results it follows that (70) man's 
E (n) into a. simnly connected domain E(P) containing the 
•point infinity and hounded hy & closed curve P . 

Prom (03) and (65) it also follows that as — £ co 

> 3(1 - p S ), 3* 0, > 0, 2Z. > 0(1 + p S ) (04) 

3 1 3r) 3 i 3rj 

i. 6 

Finally, if d£ = i dae ^ is a complex line element 
situated at a point of PI and normal to PI , the correspond- 
ing line" element dz is, hy (70) and (Cl) given hy 

dz = iO !&' [ 1 ~ 1 / n (i + p, 3 iO'[ 2 / n ) e 10 */ 11 da 

By comparison of this with (02) it is seen that the directi on 
nor mal t o PI i s taken into the dir e ct i on normal t o P . 
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M _ M + M- . cCO ' 1-1 / 31 - ll 3 &' 1 - 1 / 11 ) 


And 


dz 

3r) 


= i ("#| . ■§&)= iod?' 1 - 1 / 11 + p, a 9< l+l /“) 

v 3 & ac * 


= l m (2l Ml) = 0 a ]S'[ 3 - 2 /“ (1 - |»* |GM[ */*) 
\ 35 Sri/ 


d(x, y ) _ 

B( i i Tl) “““V 3| 3r) 



b 


4 




V 


TABLE I 

Tbe distorted speed (q*), the Mach nuBher (K) , the 
dlaeuslonleas deasity ({), and the coefficient of the symnetrized 
hodograph equation (T) as functions of the dimensionless speed 
q. (f = 1.405; = .911976; only the mantissas of log q* 

are given.) 


n 

log^q* 

' 

q* 

mm 

f 

T 

.01 

160455 

,014169 

,01000 

.99995 

1.00000 

.02 

461434 

,028956 

.02000 

.99060 

1.00000 

.05 

657472 

.043398 

.05000 

.99955 

1.00000 

.04 

762555 

.057854 

.04001 

.96920 

1.00000 

.05 

860146 

.072301 

.05001 

.99876 

1,00000 

.06 

038809 

.086750 

.06002 

.99820 


.07 

005014 

,10116 

.07005 

.98755 

.99999 

.09 

062842 

,1]M7 

.08005 

.99660 

.99999 

.00 

115609 

.12996 

.09007 

.98595 

.99998 

.10 

159550 

.14435 

.10010 

.99501 

.99997 

.11 

£00525 

.15666 

.11014 

.99396 

.99996 

.12 

£56060 

.17501 

.12018 

.00282 

.99994 

.15 

£72540 

.16750 

,13022 

.99157 

.99991 

.14 

504458 

,20158 

.14028 

.96023 

.99988 

.16 

354085 

.21582 

.15054 

.99879 

.99984 

.16 

561771 

.23002 

.16042 

.93785 

.99930 

.17 

567738 

.84420 

.17050 

.08561 

.99974 

.18 

412176 

.85635 

.18059 

.90386 

.99967 


455249 

.27243 

.19070 

.90205 

.99959 

.20 

457094 

.£8646 

.20081 

.96012 

.99950 

.21 

477650 

,30049 

.21094 

.97809 

.99959 

.22 

497557 

.51445 

.£2109 

.97507 

.99926 

.25 

516562 

.32837 

.25124 

.97376 

.99911 

.24 

534521 

.34883 

.24141 

.97145 

.99895 

.25 

551505 

.35604 

.25160 

.96904 

.99875 

.28 

667968 

.36980 

.26180 

.96654 

.P9855 

.27 

505762 

.30350 

.87808 

.96594 

.99829 

.28 

598656 

.50715 

.28225 


.99601 

.20 

613553 

.41071 

.89250 

.95647 

.99769 

.«o 

687590 

.42428 

.30277 

.95560 

.99755 


TA3IS I. f Continued) 


4 

log q* 
10 

q* 

| 


mm 

.51 

641129 

.45766 

.31306 

„ .95265 

.99695 

.32 

654211 

.45104 

.32337 

.94958 

.99652 

.53 

666833 

.46454 


.94643 

.99604 

.34 

679051 

.47758 

§v &MM 

.94519 

.99551 

.35 

690687 

.49071 


.95906 

.99492 

.36 

708845 

.50578 

.56482 

.93644 

.99427 

.37* 

713298 

.51677 

.37524 

.95294 

.99556 

,38 

784009 

.52967 

.58568 

.98934 

.99278 

.59 

734392 

.54249 

.59615 

.92566 

,99192 

.40 

744465 

.55588 

.40664 

.98189 

.99099 

.41 

754835 

.56766 

.41716 

.91804 

.90997 

.42 

763721 

.58059 

.48771 

.91410 

.96888 

.43 

778932 

.59263 

.45928 

.91003 

.98765 

.44 

781881 

.60517 

.44889 

.90587 

.96653 

.45 

790576 

.61741 

.45962 

.90178 

.96490 

.46 

799027 

.62955 

.47018 

.99750 

.98556 

.47 

807244 

.64157 

.46088 

.99515 

.96168 

.46 

815234 

.65348 

.49161 

.83872 

.97966 

.49 

825004 

.66528 

.50256 

.33420 

.97789 

.50 

650563 

.07696 

.51516 

.87961 

.97577 

.51 

857916 

.68852 

.52399 

.97494 

.97547 

.52 

345071 

.69996 

.55485 

.87019 

.97090 

.55 

852032 

.71127 

.64575 

.96557 

.96632 

.54 

858805 

.72246 

.55669 

.96047 

.96545 

.55 

865396 

.75349 

.56766 

.95549 

.96255 

.56 

871608 

.74440 

.57868 

.85044 

.95898 

.57 

E78047 

P 

.58973 

.84532 

.95558 

.58 

804116 

F (Jl 

.00065 

.84015 

.95150 

.59 

890010 


.61190 

.93406 

.94755 

.60 

896761 

.78661 

.62314 

.82953 

.94264 

.61 

901348 

.79679 

.63437 

.82412 

.95800 

.62 

906769 

.80681 

.64504 

.61865 

.95280 

,63 

912042 

.81666 

.65695 

.61313? 

,68781 

.64 

917164 

.88335 

.66852 

•B0751 

.92119 

.65 

022138 

.85587 

.67973 


.91472 

.66 

086966 

.64581 

.69119 

.79612 

.90775 

.67 

951649 

.65458 

.70270 

.79052 

.00085 

.68 

956190 

.66536 

.71426 

.78447 

.69817 

.69 . 

940590 

f «m5 

,72568 

.77855 

.38847 

•70 

. — 

944849 

.88074 

.73755 

.77258 

.67408 1 
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TABU I* (Concluded) 


q . 


q * 

Iff 

f 

T 

.71 

948969 

.08914 

.74927 

.76654 

•86396 

.72 

952951 

.89733 

.76105 

.76045 

.85303 

.73 

956795 

.90551 

.77289 

.75431 

.84121 

.74 

960501 

.91306 

.78479 

.74811 

.82844 

.75 

964070 

.92060 

.79675 

.74186 

.81459 

.76 

967500 

.92790 

.80877 

.73555 

.79957 

.77 

970791 

.93466 

.82085 

.72920 

.78325 

.78 

973942 

.94176 

.83300 

.78279 

.76547 

.79 

976951 

.94831 

.84521 

.71634 

.74606 

- .80 

979817 

.95459 

.85749 

.70984 

.72481 

.81 

982536 

.96059 

.86984 

.70329 

.70146 

.82 * 

985106 

.96629 

.88226 

.69670 

.67570 

.83 

987521 

.97168 

.89475 

.69006 

.64713 

.84 

989777 

.97674 

.90732 

.68339 

.61524 

.85 

991866 

.98145 

.91990 

.67667 

.57934 

.86 

993780 

.98578 

.93267 

.66991 

.53846 

.87 

995505 

.98970 

.94547 

.66312 

.49118 

.88 

997025 

.99317 

.95834 

.65628 

.43520 

.89 

998313 

.99612 

.97130 

.64942 

.36626 

.90 

999326 

.99845 

.98434 

.64252 

.27434 

.91 

999957 

.99990 

.99747 

.63558 

.11190 

% 

000000 

L. 00000 

1.00000 

.63425 

.00000 


TABLE II 

Constants entering in the computation of the con- 
jugate compressible flow for V — - 1. 


* 


n=iL- M 2 *, 

1/n 

/t=d*.o 



K =^r n 


.05 

1.001 

.999 

.025 

.001 

39.98 

40.16 


.10 

1.005 

.995 

.050 

.003 

19.95 

20.26 


.15 

1.011 

.989 

.075 

.006 

13.26 

13.65 

* 

• 20 

1.021 

.980 

.101 

.010 

9.90 

10.38 


.25 

1.033 

.968 * 

.187 

.016 

7.87 

8.43 

- 

.30 

1.048 

.954 

.154 

.024 

6.51 

7.13 

• 

.35 

1.068 

.937 

.181 

.033 

5.53 

6.21 


.40 

1.091 

.817 

.212 

.045 

4.72 

5.53 

* 

.45 

1.120 

.895 

.238 

.057 

4.21 

5.00 


.50 

1.155 

.866 

.266 

.072 

3.73 

4.58 


.55 

1.197 

.835 

.300 

.090 

3.34 

4.23 


.60 

1.250 

.800 

.333 

.111 

3.30 

3.95 


.65 

1.306 

.760 

.369 

.136 

2.71 

3.71 


.70 

1.400 

.714 

.408 

.167 

. 2.45 

3.50 

* 

.75 

1.535 

.651 

.454 

.206 

2.20 

3.36 

* 

.80 

1.667 

• 6Q0 

.500 

.250 

2.00 

3.17 


.85 

1.898 

.527 

.557 

.310 

1.80 

3.04 

> 

.90 



.627 

.393 

1.60 



.95 



.724 

.524 

1.38 


/ 

1.00 



1.000 

1.000 

1.00 
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Fig. 3 
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